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The electron-hole states of semiconductor quantum dots are investigated within the framework of 
empirical tight-binding descriptions for Si, as an example of an indirect gap material, and InAs 
and CdSe as examples of typical III-V and II- VI direct-gap materials. The electron-hole Coulomb 
interaction is largely insensitive to both the real-space description of the atomic basis orbitals and 
different ways of optimizing the tight-binding parameters. Tight-binding parameters that are op- 
timized to give the best effective masses significantly improve the energies of the excitonic states 
due to the better single-particle energies. However, the Coulomb interaction does not vary much 
between different parameterizations. In addition, the sensitivity of the Coulomb interaction to the 
choice of atomic orbitals decreases with increasing dot size. Quantitatively, tight-binding treatments 
of correlation effects are reliable for dots with radii larger than 15-20 A. The calculated excitonic 
gaps are in good agreement with recent photoluminescence data for Si and CdSe but agree less well 
for InAs. Further, the effective range of the electron-hole exchange interaction is investigated in 
detail. In quantum dots of the direct-gap materials InAs and CdSe, the exchange interaction can 
be long-ranged, extending over the whole dot when there is no local (onsite) orthogonality between 
the electron and hole wave functions. By contrast, for Si quantum dots the extra phase factor due 
to the indirect gap effectively limits the range to about 15 A, independent of the dot size. 

PACS numbers: 71.15.Fv, 71.24.-Fq, 71.35.-y, 71.70.Cm 



I. INTRODUCTION 

The optical and electronic properties of semi- 
conductor quantum dots have been -studied both 
experimentallyEl~llj and theoreticallyliil^tJ for a wide 
range of sizes, shapes and materials. This work is 
stimulated both by a fundamental interest in quantum- 
confined systems and by the applicability of quantum 
dots in nanoscale devices. Experimentally, significant 
recent improvements ip-|both growth techniquesc3 and 
single-dot spectroscopyn^ have enabled detailed studies of 
the energy spectra of electron-hole complexes (excitons, 
biexcitons, trions, etc.). Theoretically, the most sophis- 
ticated theocetical approaches are multi-band effective- 
mass theory,Q'll2l fpnnkjcal pseudo-potcntial,E3^L3 tight, 
binding methods ,li3~Ej and aoiasiparticle calculationsEHl 
with the GFK approximation.Eil 

Quantum dots are intermediate between molecular and 
bulk systems. This is reflected in the different theoret- 
ical approaches; effective mass theory treats the dot as 
a confined bulk system whereas pseudo-potential theory 
aims at a detailed atomistic description of the wave func- 
tions. Tight-binding theory compromises between these 
two approaches by including an atomistic description but 
limiting the local degrees of freedom to a small basis set. 



Therefore, the computationally less costly tight-binding 
method can be used to study large quantum dots, up 
to 25 nm size, without severely restricting atomic-scale 
variations in the wave functions. However, since the 
tight-binding matrix elements are empirically optimized 
without introducing any specific atomic orbitals, there 
is no direct way to calculate other matrix elements such 
as Coulomb and exchange matrix elements. Therefore, 
calculations involving electron-hole interactions require 
a selection of specific atomic orbitals that cannot be ex- 
plicitly related to the tight-binding parameters. Hence, 
a key question that we address in this work is how sensi- 
tive the calculated electron-hole properties are to specific 
choices of orbitals used to calculate the Coulomb matrix 
elements. 

We study spherical semiconductor crystallites centered 
around an anion atom in a zincblende structure. We 
choose Si, InAs, and CdSe as examples of an indirect-gap 
material, and typical III-V and II- VI direct-gap materi- 
als respectively. Experimentally, the low-lying exciton 
energjifis have been measured for dota_up to 20 Ajj'adius 
in Si^ and up to 40 A radius in InAscI and CdSe.0 

We use the— empirical nearest-neighbor sp^s* tight- 
binding modelo for the electron and hole single-particle 
wave functions. In order to calculate electron- hole 
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Coulomb and exchange matrix elements, we describe the 
real-space atomic basis orbitals s,Px,Py,Pz, and s* with 
Slater orbitals as a starting point. Both the Coulomb and 
exchange interaction are screened by a dielectric function 
depending on both dot size and the distance between the 
particles. The energies of the electron- hole states are 
obtained by diagonalizing the configuration interaction 
matrix generated by the lowest-lying electron and hole 
states with typical convergence of a few meV. We exam- 
ine the sensitivity of the electron-hole energies on both 
the choices of atomic basis orbitals and the tight-binding 
parameters. 

Within the tight-binding model, the single-particle 
Hamiltonian can be improved by either increasing the 
number of basis orbitals or including interactions between 
more distant-atoms. An s* orbital was first introduced 
by Vogl et a/.Lj to improve the conduction bands near the 
X point. To some degree, an s* orbital can mimic band 
structure effects that should be attributed to d bands. 
However, it is Itely that d orbitals need to be added in 
some materialsO An_alternative is to add next-nearest 
neighbor interactionai3 within the sp^s* model, which 
could improve the band structure without increasing the 
computational cost of generating Coulomb matrix ele- 
ments, since the number of orbitals is unchanged. Never- 
theless, for the topics we discuss in this work the nearest- 
neighbor sp'^s* model is a good starting point. We here 
focus on how the single-particle energies can be improved 
by tight-binding parameters specifically optimized to give 
good effective masses, and on how single-particle wave 
functions from different parameter sets affect the two- 
particle Coulomb interactions. 

One interesting issue concerning quantum dots is 
the range of thcr-electron-hole exchange interaction. 
Franceschetti et alSB show that due to the lack of local or- 
thogonality between electron and hole wave functions the 
exchange interaction can extend over the whole dot. We 
investigate in detail the effective range of the exchange 
interaction by applying a cutoff range to the Coulomb 
potential. The origin of the characteristic range of the 
exchange interaction is revealed by the analysis of the 
'exchange charge density' of the electron-hole pair. 



II. THEORY 

A. Hamiltonian of an electron-hole pair 

The effective Hamiltonian of an electron-hole pair con- 
tains a single-particle term and an electron-hole term as 
follows (excluding spin-orbit coupling): 

= ^ I dV ^+(r) {-iv2 + y(r)} ^,(r) 

+ -E //d'^' ^^(r')V3+(r)V3.(r)Vv(r') 



where 
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(2) 



with the spin index s and the tight-binding orbital index 
n, which includes atom-site index i and orbital-type index 
7. The functions 0n(r) are tight-binding basis orbitals, 
which in this work are s,Px,Py,Pz or s*. The operators 
and h'X s are annihilation and creation op- 
erators of an electron and a hole for the tight-binding 
basis orbital (/)n(r), respectively. Note that a„^s and 
are applied to conduction electron states, while 6„^s and 
6+^ are applied to valence electron states. The dielectric 
function e(|r' — r|, i?) is assumed to be a function of both 
the separation |r' — r| of two particles and the dot radius 
R. Note that we use atomic units for all the equations in 
this paper. 

The single-particle potential V is implicitly defined 
through the empirical tight-binding matrix elements. We 
use the nearest-neighbor sp^s* tight-binding description. 
The structure of a quantum dot, is modeled as an anion- 
centered zincblende structure.E3 Dangling bonds on the 
surface are removed by explicitly shifting the energies of 
the corresponding hybrids well above the highest calcu- 
lated electron states. This treatment mimics a dot whose 
surface is efficiently passivated with, for example, hydro- 
gen or ligand molecules. 

By multiplying the electron and hole eigenstates from 
the solution of the tight-binding Hamiltonian, we obtain 
an electron-hole basis set \eh)\js,ms) where \js,rns) is 
an electron-hole spin state, i.e., either the singlet compo- 
nent |0,0} or one of the triplet components |1,0), 
or |1, — 1). The spatial part \eh) is the product of an elec- 
tron and a hole wave function that in real space has the 
form: 

V'e(re)V',t(r/i) = Ce,mCln^<pmire)(Pn^{rh)- (3) 



This description closely follows that of Leung and 
Whaley.B 

The single-particle Hamiltonian can be written in 
terms of the electron-hole basis set and its eigenvalues: 

J^singio = ^ {Ee - Eh)\eh)\jsms)(jsms\{eh\, (4) 

ehjs rris 

where and Eh are the electron and hole energies of 
the single-particle Hamiltonian. 

Projecting the electron-hole Hamiltonian onto the two- 
particle basis set yields the electron-hole Hamiltonian 
with a CM)£>pib interaction J and an exchange inter- 
action Kxl 



e(|r' - r|,i?) |r' - r| 



, (1) 



(5) 
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e'h' eh 



e(|r'-r|,i?) |r'-r| ' 

K = 2Sj^ \e'h'){eh\ 



e' h' eh 



e(|r'~r|,i?) |r'-r| ' 



(6) 



(7) 



where is unity for a singlet state and zero for a triplet 
state. The factor 2 in front of in Eq. ^ is due to 
the fact that the exchange interaction allows two final 
electron-hole spin states | te, ih) and | J.^, f/j) for an ini- 
tial state I te, ih) (or | |e, T/i))- In contrast, the Coulomb 
interaction requires that the spin of the electron should 
be the same between an initial and a final state, as should 
the spin of the hole. 

The Coulomb interaction J describes the scattering 
of the electron from e to e' and the hole from h to h' , 
whereas the exchange interaction K describes the recom- 
bination of a pair e,h at r and the recreation of a pair e',/i' 
at r'. Since we do not include spin-orbit couplings, the 
spin states are not coupled to one another in the present 
Hamiltonian. Therefore, we will use only the spatial part 
\eh) of the electron-hole basis set from now on. The only 
constraint that the spin state gives to the Hamiltonian is 
the spin-selection rule in K . 

The matrix elements of the electron-hole interaction 
Hamiltonian can now be rewritten in terms of integrals 
over the tight-binding basis orbitals by replacing ?/'e(re) 
■with Z]„Ce;„(^„(re) and -0/i(rh) with X]„ c/i;„0„(r,i): 

{e'h'\J\eh) 

= - X! Ce';„iCe;„2 4.„3C?i/;„4 u;(ni, 722; ^3, rii), (8) 

{«} 
{e'h'\K\eh) 

uj{ni,n2;n3,n4), (9) 

{n} 



where 



w(ni,n2;"3,«4) 



>;^(r')</.„,(r')C3(r)</'n4(r) 



(10) 



e(|r' -r|,i?) |r' - r| 

Following Leung and Whaleyjl^ we approximate the 
Coulomb and exchange interaction matrix elements by 
considering only terms having at most two distinct ba- 
sis orbitals. This approximation is reasonable since the 
integrals involving more than two different odjjtals are 
typically small compared to the kept integrals .E3'llZl Inte- 
grals in Eq. (^) with ni = and 713 = 714 are Coulomb 
integrals and those with ni — 77,4 and n2 = 713 or with 
ni = 77.3 and 712 — are exchange integrals. Hence, we 
define: 



wcoui(ni,7i2) = uj{ni,ni;n2,n2) 
t^cxch("-i, f^2) = uj{ni,n2;n2,ni) 
= cj(r2i, 772; 711,712), 



(11) 
(12) 



where the equality of the two exchange integrals is valid 
when the tight-binding orbitals are real, as they are in 
this work. 

To make our notations clear, note that the Coulomb 
and the exchange integrals of the basis orbitals are the 
interactions between the tight-binding basis orbitals. By 
contrast, the Coulomb and the exchange interactions J 
and K are interactions between the electron and hole 
wave functions. In fact, the Coulomb interaction has 
contributions from both the Coulomb and exchange in- 
tegrals as does the exchange interaction. 

Finally, we can describe the electron-hole matrix ele- 
ments in terms of the Coulomb and exchange integrals. 

{e'h'\J\eh) 

= - X! '^e';niCe;ni4;„2C/i';„2 ^Coul (fll , "-2 ) 
ni ,n2 

~ C*/.„jCe;Ti2C/^;TiiC/i';r!2 ^cxch (^^ 1 , '^■2 ) 

~ 51 C*,.„^Ce;„2C;;.„^CV;„i t^cxch(n-l,ri2), (13) 
Tii,n2 

{e'h'\K\eh) 

= '^^js X! '^e';niC'i';«iC;;.„^Ce;„2 ^Coul (»^1 , »^2) 
:ni 

ni ,n,2 

+ E Ce';„iCft';n24;«iCe;n2 ^oxch ("•! , "2 ) • (14) 

"1 ,"2 

To evaluate these integrals, we need a real-space de- 
scription of the tightJpinding basis orbitals. To start, we 
follow Martin et a/£3 and model j-the tight-binding or- 
bitals with atomic Slater orbitals. E3 The Slater orbitals 
are single-exponenlial functions with the exponent given 
by the Slater rulesEZi designed to yield a good approxi- 
mation of the effective radius and the ionization energy. 
Especially, an s* orbital is modeled as an excited s orbital 
by promoting one valence electron to the s orbital of the 
next shell. The Slater orbitals are an arbitrary choice in 
the sense that they are not explicitly related to the tight- 
binding parameters. However, as shown in Section IH, 
even in small dots the electron-hole Coulomb interaction 
is not very sensitive to variations in the orbital integrals. 

Onsite Coulomb and exchange integrals, in which both 
orbitals are centered on the same atom, are calculated us- 
ing a Monte Carlo method with importance sampling for 
the radial integrations. The uncertainty of the Monte 
Carlo results is within 1%. The angular part is treated 
exactly by expansion in spherical harmonics. However, 
the Appendix shows that the integral values must be con- 
sidered to be uncertain to about 20-30% due to the ar- 
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bitrariness of the orbital choice and the effects of orthog- 
onaUzation. 

Offsite exchange integrals, where the two orbitals are 
centered on two different atom sites, are negligible even 
for nearest-neighbor integrals. These offsite exchange in- 
tegrals decrease quickly as the distance between atom 
sites increases, due to the localization and orthogonality 
of the orbitals. In particular, we show in the Appendix 
that even nearest-neighbor exchange integrals are neg- 
ligible as an effect of orthogonalization between offsite 
hybrids. 

Regarding offsite Coulomb integrals, Leung and 



Whaley— estimate these integrals with 
formulaEa modified to include screening: 

t^Coul(", n) EE UJCon\[il, i'l) 

1 1 



the Ohno 



F, (15) 



where and Rj' are atom site vectors. ^^^^^[(^7, 17') is 
an unscreened onsite Coulomb integral. The superscript 
designates an unscreened quantity. For the case of bi- 
nary compounds, '^coui(*7' replaced by the average 
over the two different atoms. The integrals are screened 
by the dielectric constant e(|Ri — Ri'|,i?). This screen- 
ing is the only modification to the original Ohno formula. 
From here on, we will refer to this modified Ohno formula 
in Eq. simply as the Ohno formula. 

To test the validity of the Ohno formula in the case that 
two orbitals are on close atom sites, we calculated tke, 
offsite Coulomb integrals with a Monte Carlo methodo 
within 1% uncertainty and compared these values with 
those from the Ohno formula. The Ohno formula severely 
underestimates the offsite integrals as the distance be- 
tween two atom sites becomes small. For example, the 
Coulomb integral between the bonding orbitals (see be- 
low) from nearest neighbors in a Si quantum dot with 
radius 18.9 A given by the Ohno formula is 0.58 eV, 
while the Monte Carlo calculation gives 2.35 eV. For 
next-nearest neighbors, the Ohno formula and the Monte 
Carlo calculation give 0.38 eV and 0.58 eV, respectively, 
and for the third-nearest neighbors 0.33 eV and 0.36 eV. 

The reason that we obtain a big discrepancy between 
the Ohno formula and the Monte Carlo calculation for 
the bonding-orbital integrals is that the effective distance 
between the bonding-orbitals is smaller than the spacing 
between the nearest-neighbor atom sites. In fact, the 
spatial overlap of the bonding-orbitals is as big as that 
of the orbitals on the same atom site. In addition, the 
spatial dependence of the dielectric function is not fully 
taken into account in the Ohno formula. This approxi- 
mation becomes critical when the range of variations in 
the dielectric function is comparable to the effective dis- 
tance between orbitals. In that case, the effective dielec- 
tric function cannot be represented by e(|Ri — Rj|,i?). 
Therefore, we use the Monte Carlo values for the onsite 
and the nearest-neighbor integrals and the Ohno formula 
for the rest of offsite integrals. For clarity, we summarize 



the methods for the computation of the Coulomb and 
exchange integrals in Table |. 

We use a size and distance dependent dielectric func- 
tion to screen the Coulomb and exchange interaction of 
the electron-hole pair. The dielectric function, as a func- 
tion of the separation r of two particles and the radius R 
of a quantum doJ— is approximated by the Thomas- Fermi 
model of Restatl-and the Penn model generalized for 
quantum dots.Ell'Ej The separation dependence is given 
by the Thomas-Fermi model, while the size dependence 
is given by the Penn model. This way of combining the 
two models to describe the dielectric function is taken 
from Rcf. 13 



e(r,i?) 



where 



e^*(i?) gro/[sinhq(ro - r) + qr] , r < tq 



Aot 



(R) 



hulk 



1) 



f nibulk 



A)2 



(16) 



(17) 



The Thomas-Fermi wave vector q is (4/7r)^/^(37r^no)^/^, 
where the valence electron density no — 32/ao in a 
zincblende structure. The screening radius ro is deter- 



mined by the condition sinh grg / qr^ 



^dot 



(R). The shift 



A = i?2 — -E-gap , where E2 is the energy of the first pro- 
nounced peak in the bulk absorption spectrum. The en- 
ergies -Egap'^ and E^°p{R) are the single-particle gaps for 
bulk and a dot with radius R, respectively, ej^"^ is the 
dielectric constant for the bulk material. 

The unscreened onsite Coulomb and exchange inte- 
grals for the sp^s* basis set are listed in Table ||, and 
the screened onsite Coulomb integrals are listed in Ta- 
ble [II. The screened nearest-neighbor Coulomb inte- 



grals are listed in Table IV. The integrals with s and p 



orbitals are calculated using the four hybridized orbitals 
along bonding directions, defined as: 

\spl)^^{\s)±\p,)±\Py)±\p,)), 

\spl) = ^(|s) ± \Px) T \Pv) T \Pz)), 

\spI) ^^{\s)t\Px)±\Pv)t\Pz)), 

\spl)^li\s)T\p.)T\Pv)±\Pz)), (18) 

where the upper sign is for an anion and the lower sign 
is for a cation. 

As shown in Tables |l| and [II , the screening effect is 
significant even for onsite integrals. This is because the 
screening radius ro is 2-4 A and is similar to the ef- 
fective radius of the tight- binding basis orbitals. The 
comparison of the unscreened and screened onsite inte- 
grals shows that the effective screening of these integrals 
is generally about half the long range screening given by 
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e^*(i?). Further, based on this observation, we use half 
of the long-range dielectric constant to screen the onsite 
exchange integrals that we-did not calculate explicitly by 
the Monte Carlo mcthodO 



B. Lowest excitonic states 

To obtain the excitonic states near the band edge, we 
diagonalize the configuration-interaction matrix in the 
\eh) basis set given by the sum of i?singic and i?c-ii de- 
fined by Eqs. (Q)- (^. We include sufficient electron and 
hole states in the configurations to achieve convergence 
of the first few excitonic states to within a few meV. 
The typical number of electron and hole states needed is 
about f0-f5 each. 

There are two types of Hamiltonians depending on the 
total spin of the electron-hole pair. The Hamiltonian for 
a spin singlet includes both the Coulomb and the ex- 
change interaction, whereas the Hamiltonian for a spin 
triplet has only the Coulomb interaction. By diagonal- 
izing these two Hamiltonians separately, we obtain a set 
of spin-singlet and spin-triplet excitonic states. The low- 
est excitonic state is the lowest triplet state due to the 
absence of the positive exchange interaction. However, 
since only spin-singlet states are optically allowed the op- 
tical excitonic gap is the energy of the lowest spin-singlet 
state. 

This order of the states can be seen most easily by 
applying first-order perturbation theory to the electron- 
hole pair made from the highest hole state and the lowest 
electron state, which yields 

^^singiot ^E,-Eh + {eh\J\eh) + {eh\K\eh). (19) 
Et,ipi,t^E,-Eh + {eh\J\eh). (20) 

The sign of {eh\J\eh) is always negative and the sign 
of {eh\K\eh) is always positive. Therefore, the en- 
ergy of the lowest spin-triplet excitonic state is smaller 
than the energy of the lowest spin-singlet excitonic 
state by {eh\K\eh) within first-order perturbation the- 
ory. The energy difference between these excitonic states 
is the exchange splitting. Further, we denote the dif- 
ference between the lowest spin-triplet energy and the 
single-particle energy gap as the Coulomb shift, which 
is {eh\J\eh) in first-order perturbation theory. The 
Coulomb shift is the main correction to the single-particle 
gap since the Coulomb interaction is roughly one order 
of magnitude larger than the exchange interaction. 

This simple description becomes more complicated 
when configuration interaction is included due to the cor- 
relation of several electron-hole configurations near the 
band edges. However, the main ideas of Coulomb shift 
and exchange splitting are still valid. Generally, as we 
include more electron-hole configurations the Coulomb 
shift and the exchange splitting increase and converge. 



C. Effective range of the exchange interaction 

The long-range component of the exchange interac- 
tion was investigated by Franceschetti et alBB with the 
pseudo-potential method. They show that there is a long- 
range component in the monopole-monopole exchange 
interaction for several direct gap semiconductor quan- 
tum dots. To verify this long-range exchange interaction 
with the tight-binding model, we follow their approach 
using a cutoff potential. With the step function 0(r), 
we replace the Coulomb potential with a cutoff poten- 
tial 8(rc — |r' — r|)/|r' — r| to obtain an exchange inter- 
action that is a function of the cutoff distance Tc. 

The unscreened exchange interaction with the cutoff 
potential for electron- hole state \eh) is 

{eh\K"{rc)\eh) 

^ J I ^3,, ,3,C(rO^.WJ(r)^e(r) ^^^^ _ 

cl.ni^h;niCl^-n2^e-,ri2 ^Coul(^l'^2) 

rii .712 

X e(rc - |R„i - R„J) 

+ 4 ^ C*.„iCh;n24;n2Ce;«i W°^eh ('^l . "2 ) ■ (21) 

"1 ,"2 

where the superscript refers to the unscreened inter- 
action. In line with the discrete spatial character of the 
tight-binding model, we make an approximation that re- 
places the true cutoff potential with one based on the site 
indices Q{rc — |R,ii — RnaD- If there is a long-range ex- 
change interaction, it would stem from the first term of 
Eq. ( ^l| ) which includes the Coulomb integrals. The sec- 
ond term, the sum of exchange integrals, has only onsite 
integrals, since all offsite exchange integrals are negligible 
as shown in the Appendix. 

To understand the physical origin of the long-range 
exchange, we expand the exchange interaction i^T" in a 
multipole expansion: 

{eh\Krh) « E '^^ii'^'nt' + " ^^'"^ 

-I- onsite interaction. (22) 

Here the 'exchange charge density' q{'Ri)eh at atom site 
Ri is a monopole moment defined as 

= 5]ce;^(R.)c;;;y(R.) I d'r ct>^^{r^*y{r) 

= Ce;7(Ri)c^.^/(Rj)^77' 

77' 

= Eq;^(r,)c;;.^(r,), (23) 

7 
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where J dfii is defined to integrate only the orbitals on 
the atom site Ri. For clarity, the tight-binding orbital 
index n is specifically replaced by the atom-site index i 
(or Ri) and the orbital- type index 7. Note that the final 
expression for q{'Ri)eh has a sum over only one orbital- 
type index due to the assumed orthogonality of the tight- 
binding basis orbitals. The distribution of the exchange 
charge density determines the long-range character of the 
exchange interaction. If the exchange charge density is 
zero, that is, the electron and hole states are locally (on- 
site) orthogonal, there is no exchange interaction beyond 
the onsite contribution according to Eq. (p2|). In con- 
trast, if the exchange charge density is nonzero due to 
the local nonorthogonality of the electron and hole wave 
functions from site to site, a long-range exchange inter- 
action is caused by the monopole-monopole interaction. 



III. RESULTS 

A. Real-space description of basis orbitals 

The empirical tight-binding model has an inherent dif- 
ficulty concerning the tight-binding basis orbitals. The 
real-space description of the basis orbitals is not provided 
since the tight-binding matrix elements are determined 
by fitting to the bulk band structure. However, to include 
electron-hole correlations the electron-hole Coulomb and 
exchange matrix elements need to be calculated, which 
requires an explicit choice of real-space basis orbitals. 
Since this choice is largely arbitrary in the sense that 
there is no way to connect the chosen basis orbitals to 
the empirically chosen tight-binding parameters, we need 
to test to what degree the choice of orbitals affect the 
electron-hole Coulomb interaction. 

We perform this test by scaling the onsite Coulomb 
and exchange integrals from the values listed in Tables || 
and III. This scaling scheme is an indirect way of test- 
New 



ing the sensitivity on the real-space description, 
offsite Coulomb integrals are calculated by replacing the 
unscreened onsite integrals with the scaled ones in the 
Ohno formula, Eq. (p^. Note that the offsite Coulomb 
integrals are not directly scaled by the same factor as the 
onsite integrals, but change only indirectly through the 
scaled unscreened onsite integrals in the Ohno formula. 
Therefore, this scheme is in effect changing atomistic de- 
tails of the basis orbitals. 

Figure |l| shows the variation of the Coulomb interac- 
tion between the highest hole state and the lowest elec- 
tron state with the scaling of the onsite Coulomb and ex- 
change integrals. It shows that as the dot size increases, 
the sensitivity of the Coulomb interaction to the onsite 
integrals decreases. For example, if the onsite integrals 
are reduced by 50%, the reduction in the Coulomb en- 
ergy varies from 20% in the smallest shown dot to only 
5% at 30 A radius. Since the contribution from the onsite 
integrals decreases as the dot size increases, the specific 



model of the real-space functions for the basis orbitals is 
less critical for larger dots. 

We can explain this effect by a closer look at the Ohno 
formula for the offsite integrals. In the limit of large 
distances between two atom sites, the offsite integrals 
become independent of the onsite integral values. The 
unscreened offsite integral in the limit of large distance 
between the two atom sites is 



'^Coul(«7,*V) 

1 

^ |R,:-R,:'I 



2wcoui(*7,*y)' IR.-R. 



(24) 



Therefore, the offsite integrals become a point-charge in- 
teraction, making the atomic-scale details of the basis 
orbitals irrelevant in this limit. 

The significance of Figure |] is that it quantifies this 
qualitative explanation for decreasing sensitivity with in- 
creasing dot sizes. For example, with a targeted 10% ac- 
curacy in the Coulomb interaction, only a 50% accuracy 
for the onsite integrals is needed for dots of radius larger 
than 30 A, while for dots of 10 A radius only 20% error 
can be afforded in the onsite integrals. The discussion in 
the Appendix shows that the dominant integrals must be 
considered uncertain to about 20-30%. The tight-binding 
description of correlation effects can therefore be consid- 
ered reliable for dots with radii larger than 15-20 A. 



B. Excitonic states near the band edge 

We apply the tight-binding configuration interaction 
scheme described in Section 



quantum dots in order to ca 



11 A to Si, InAs, and CdSe 



culate the lowest excitonic 



states near the band edge. We have included sufficient 
electron and hole states to converge the energies to within 
a few meV. 

For Si, pWe first used the tight-binding parameters of 
Vogl et aZo to determine the single-particle Hamiltonian 
matrix elements. As shown in Figure]^, the excitonic gap 
using their tight-binding parameters gives a discrepance 
as large as 0.3 eV compared with experimental data.cl 
The tight-binding parameters of Vogl et al. necessarily 
give poor effective masses since the parameters are de- 
termined by fitting to the energies of only the F and X 
points in the bulk band structure. The resulting effective 
masses are listed in Table 0. The comparison with the 
experimental values listed in Table ^ shows that their 
parameters fail to produce good effective masses espe- 
cially for the transverse effective masses at the lowest 
conduction energy near X . 

To improve the effective masses, we replace the pa- 
rameter set of Vogl et al. with our parameter set listed 
in Table VI. Our parameter set is optimized with a ge- 
netic algorithm by fitting the effective masses as well as 
the energies at high symmetry points of the bulk band 



structurcta The resulting effective masses are listed in 
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Table 0. One important note is that we use two dif- 
ferent parameter sets to separately optimize the electron 
and hole single-particle states. Good effective masses are 
impossible to obtain simultaneously for both the con- 
duction and the valence band of Si with one set of pa- 
rameters within the nearest-neighbor sp^s* tight-binding 
model (see Ref. |8|). Consequently, the electron and hole 
single-particle wave functions, being generated from dif- 
ferent Hamiltonians, are not orthogonal. However, even 
though the orthogonality has not been enforced, the over- 
laps between the different electron and hole wave func- 
tions are at most 0.001. Thus, we can use these two 
different parameter sets to verify how important a role 
the effective masses play in the electronic properties of 
the quantum dots. 

Figure ^ shows the improved excitonic gaps with our 
parameter set. To further examine the effect of changing 
parameters, we can compare the electron-hole interaction 
energies with our parameters to those energies obtained 
with the parameters of Vogl et al. In particular, we com- 
pare the Coulomb shift, the energy difference between 
the single-particle gap and the lowest triplet excitonic 
energy. As shown in the inset in Figure ||, the Coulomb 
shifts from the two parameter sets are very similar. This 
insensitivity indicates that the better description of the 
excitonic gap with our parameter set is mainly due to the 
better single-particle eigenvalues and not from a change 
in the Coulomb matrix elements. 

To study direct gap semiconductors, we choose InAs 
and CdSe spherical quantum dots. The InAs tight- 
binding parameters are generated using the genetic algo- 
rithm approach,tj fitting band gaps and effective masses 
at r as well as possible, but neglecting sp in-o rbit cou- 
pling. These parameters are listed in Table VII. The re- 
sulting effective masses with these parameters are nic — 
0.024, nivi = 0.025, and m^h — 0.405, where rric is the ef- 
fective mass of the lowest conduction band at F, and rriyi 
and rjivh are defined as in Table 0. The tight-binding 
parameters for CdSe are taken from Ref. 

Figures ^ and ^ show the resulting excitonic gaps ver- 
sus the dot radius. For CdSe quantum dots, our ex- 
citonic gaps are in good agreement with^nptical gaps 
measured by photoluminescence excitationO (PLE). We 
also plot the energy gaps measured by scanning turt^ 
neling spectroscopy (STM) on a single quantum dot.EI 
The STM gaps are obtained from the difference between 
the first prominent peaks of the tunneling dl/dV spec- 
tra with positive and negative bias voltages, respectively. 
Since the STM experiment applies bias voltages to add 
or subtract electrons from the quantum dots, this ex- 
periment measures quasiparticle energies. For a finite 
system, these quasiparticle energies include the (posi- 
tive) polarization energy between the particle and the 
image charges on the surface. The polarization energy 
is roughly 2(l/eout — l/cin)/^- The smaU difference be- 
tween our single-particle gaps and the STM gaps there- 
fore suggests that the dielectric constant Cout of the sur- 
rounding material is relatively close to the dielectric con- 



stant eiii in-the dots. The results of pseudo-potential 
calculational3 are also plotted in the figures for compar- 
ison. 

For InAs quantum dots, there is a significant discrep- 
ancy as large as 0.2 eV between our sp^s* tight-binding 
excitonic ga^s and PLE gaps.Q Eight-band effective-mass 
calculationsQ and pseudo-potential calculationsEdl also fail 
to describe the experimental data, especially the lack of 
significant curvature. The recent results of Allan et al. 
show that the inclusion of d orbitals and spin-orbit inter- 
action raises the sp^s* results by almost the needed 0.2 
eV. However, their results do not include the Coulomb 
shift and should therefore be shifted down by 200-50 meV 
as the dot size increases. Recent STM measurements are 
also plotted in Fig. ^. It is consistent with the larger 
dielectric constant of InAs that the STM results in this 
case is well above the other curves by an amount similar 
to the Coulomb shift. 



C. Effective range of the exchange interaction 

One of the interesting issues related to the exchange 
interaction of the electron-hole pair is its effective range. 
Motivated by the work of Franceschetti et aZ.j2a we cal- 
culate the unscreened exchange interaction defined in 
Eq. (|2^) as a function of the cutoff distance rc to de- 
termine the effective range of the exchange interaction. 
As the cutoff distance increases for a given electron-hole 
configuration, the exchange interaction eventually satu- 
rates to a final value. If this saturation occurs over just 
a few atomic sites, we call it short-ranged, while long- 
ranged exchange implies that the saturation occurs over 
distances comparable to the dot size. For Si, Figure ^ 
shows that for the configuration with the highest hole 
state and the lowest electron state there is a region of 
strong oscillations below a cutoff distance of 15 A. The 
strong oscillations are due to the phase difference be- 
tween the electron and hole states stemming from their 
different locations in fc-space for an indirect-gap material. 
The oscillations die out beyond a cutoff distance of about 
15 A, suggesting that the effective range of the exchange 
interaction in Si quantum dots is around 15 A regardless 
of dot size. This short-ranged and oscillatory behavior is 
universal within the configurations near the band edges. 

For the direct-gap InAs and CdSe quantum dots, we 
calculate the unscreened exchange interaction for several 
of the lowest electron-hole configurations. We label the 
electron and hole states by the dominant angular mo- 
mentum character of their 'envelope functions.' Here, 
the envelope function is defined to be the coefficient of 
the dominant basis orbital. The s and p basis orbitals are 
typically dominant in the electron and hole states, respec- 
tively. In our calculations, a p-like holecj is the highest 
hole state and an s-like hole is the second-highest hole 
state. Ihis order is opposite to that of pseudo-potential 
theory.Ej However, it is possible that the spin-orbit cou- 



7 



pling, which is not included in this work, can affect the 
order of these hole states. 

As shown in Figures || and ^ direct-gap quantum dots 
show a qualitatively different behavior of the exchange 
interaction with respect to the cutoff distance from the 
behavior for Si. First, since there is no overall phase dif- 
ference between the electron and hole states, there is no 
region with oscillations for small cutoff distances. Sec- 
ond, the exchange interaction for a particular electron- 
hole pair can grow continuously up to the dot radius. 
The figures show that the exchange interaction of direct- 
gap materials is generally long ranged, extending over 
the whole dot. 

To understand why some electron-hole configurations 
have a slowly varying long-range exchange interaction, 
we analyze the long-range cornponent by a multipole- 
expansion as written in Eq. (|22|). The leading term 
of the long-range exchange interaction is the monopole- 
monopole interaction. Therefore, the distribution of the 
monopole moment, or the 'exchange charge density' de- 
fined in Eq. (23), determines the range of the long-range 
exchange interaction. 

The exchange charge density of an electron-hole pair 
has zero total charge due to the orthogonality between 
the electron and hole wave functions. There are two ways 
to satisfy this condition: the electron and hole states 
are either locally orthogonal, which is enforced in effec- 
tive mass theory due to the orthogonality between the 
Bloch functions of the valence and conduction bands, 
or globally orthogonal which is possible in the atomistic 
pseudo-potential and tight-binding theories. If the for- 
mer is true, the exchange charge density would be zero 
at each site and there would be no monopole-monopole 
interaction. That would make the exchange interaction 
of the electron-hole pair short ranged. By contrast, with- 
out onsite orthogonality the exchange charge density has 
nonzero values, causing monopole-monopole interactions 
that lead to significant long-range exchange interactions. 

To show that the character of the orthogonality of the 
electron-hole configuration determines the long-range be- 
havior of the exchange interaction, we plot in Figure ^ 
the exchange charge density of two electron-hole config- 
urations in CdSe that have a long-range exchange inter- 
action in Figure ^. Figure ^ shows the exchange charge 
density of (a) the s-like electron and s-like hole config- 
uration, and of (b) the p-like electron and s-like hole 
configuration in a plane going through the center of the 
dot for CdSe with radius 20.1 A. This fi gure shows that 
there is no local orthogonality between the electron and 
hole wave functions. The orthogonality of the electron 
and hole wave functions are instead satisfied by a p-like 
global oscillation (case a) or a 2s-like global oscillation 
(case b). These shapes of the global oscillations explain 
why the exchange interaction has growing and decaying 
regions over global distances as shown in Figures I and 0. 
Those electron-hole configurations that do not have a 
long-range exchange interaction have a much smaller ex- 
change charge density than those configurations that do 



have the long-range exchange interaction. These results 
show that local nonorthogonality of the electron and hole 
wave functions leads to a strong monopole-monopole in- 
teraction, and that the global variations in the exchange 
interaction depend on the particular way in which the ex- 
change charge density globally sums to zero for a specific 
electron-hole configuration. 



IV. SUMMARY 

We use tight-binding wave functions to calculate 
electron-hole states near the band edge for both direct- 
gap and indirect-gap quantum dots. First, we examined 
to what degree the model of the real-space atomic ba- 
sis orbitals affects the electron-hole Coulomb interaction. 
We find that the sensitivity of the Coulomb interaction to 
the real-space description of the basis orbitals decreases 
quickly as the dot size increases. Our results shows that 
tight-binding descriptions of electron-hole Coulomb in- 
teractions in quantum dots should be reliable for dots 
larger than about 15-20 A radius even for simple models 
for the basis orbitals. More detailed calculations of basis 
orbitals are required for smaller dots. 

For excitonic gaps, we obtained good agreement with 
recent experiments for both Si and CdSe quantum dots. 
However, the gaps for InAs quantum dots agree less well 
with experiment. Especially for Si, we improved the 
agreement with experimental data for the excitonic gap 
by optimizing the tight-binding parameters to give better 
effective masses compared to the parameters of Vogl et 
al. We also showed that, in contrast to the electron and 
hole single-particle energies, the electron-hole Coulomb 
interaction is not very sensitive to the choice of parame- 
ters. 

Finally, we studied the effective range of the exchange 
interaction. Replacing the Coulomb potential with a cut- 
off potential, we demonstrated the dependence of the ex- 
change interaction on the cutoff radius. The existence 
of a global exchange charge density of an electron-hole 
pair, due to the lack of onsite orthogonality, causes the 
exchange interaction to be long ranged in direct-gap ma- 
terials. For indirect materials, the calculations for Si 
show that the exchange interaction is oscillatory and has 
a range of about 15 A. 
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APPENDIX: 

Since we use empirical tight-binding wave functions, 
the choice of specific atomistic orbitals for matrix- 
element calculations is largely arbitrary. The results pre- 
sented in the main text were based on Coulomb and ex- 
change integrals calculated 'V¥ith Slater's atomic orbitals, 
obtained from Slater's rules. cll In addition, we neglected 
all offsite exchange integrals. An alternative, for un- 
screened integrals, is to use one of the standard Quantum 
Chemistry Gaussian-based commercial packages. How- 
ever, screened matrix elements cannot be obtained in this 
way, since there is no way in these codes to include a spa- 
tially varying screening function. 

Although we cannot obtain screened integrals, two im- 
portant questions can be answered by a comparison be- 
tween integrals from Gaussian- type orbitals (GTO) and 
Slater's orbitals (SO): what is the typical variation in 
the integral values for two reasonable choices of orbitals; 
and what is the effect of using nonorthogonal bond hy- 
brids rather than properly orthogonalized hybrids? The 
underlying assumption in the tight-binding approach is 
that the orbit als on different sites are orthogonal. 

Tabic |VII]| shows a comparison between orthogonal 
GTO (0-GTO), nonorthogonal GTO (NO-GTO), and 
nonorthogonal SO (NO-SO) for onsite Coulomb and ex- 
change integrals. Typically, the NO-SO and NO-GTO 
Coulomb integrals differ by 10%, whereas the (an or- 
der of magnitude smaller) exchange integrals differ by 
20-50%. Orthogonalization generally gives an additional 
10% change. The use of nonorthogonal Slater orbitals 
can therefore be estimated to imply 20% overall uncer- 
tainty in the onsite integrals. 

A similar com par ison for nearest-neighbor integrals is 



shown in Table IX. Here the difference between NO- 
GTO and NO-SO is less than 10%, but orthogonaHza- 
tion can yield a lowering of up to 30% in the Coulomb 
integrals between bonding orbitals. The most dramatic 
effect, however, is that the exchange integrals essentially 
become negligible when orthogonalized hybrids are used. 
Notably, nonorthogonal hybrids cannot be used for the 
bonding-bonding offsite exchange integrals, since these 
integrals are quite large without orthogonalization but 
are reduced by a factor of 20-30 after orthogonalization. 

In conclusion, the dominant Coulomb integrals ob- 
tained from the Slater orbitals can be considered accurate 
only to 20-30% due to the sensitivity to different func- 
tional representations and to effects of orthogonalization. 
Further, proper orthogonalization reveals that all offsite 
exchange integrals can be neglected, including those be- 
tween bonding hybrids. 
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TABLE I. Methods for the computations of the Coulomb 
and exchange integrals with respect to a site-to-site distance. 
NN stands for nearest neighbors. 





onsite 


NN 


beyond NN 


Coulomb 


Monte Carlo 


Monte Carlo 


Ohno formula 


exchange 


Monte Carlo 


neglected 


neglected 



TABLE IL Onsite unscreened Coulomb integrals, Eq. (Pt|), 
and exchange integrals, Eq. (jl^, of the sp^s* basis set, 
in units of eV. Integrals for the sp^ orbitals are calculated 
based on the hybridized orbitals along the bonding directions, 
Eq. 

Element Integral {spl,spl) [spl,spl) {spl,s*) {s*,s*) 



Si 
Si 
In 
In 
As 
As 
Cd 
Cd 
Se 
Se 



Coulomb 
exchange 
Coulomb 
exchange 
Coulomb 
exchange 
Coulomb 
exchange 
Coulomb 
exchange 



11.88 
11.88 
7.82 
7.82 
12.13 
12.13 
6.59 
6.59 
12.85 
12.85 



8.49 
0.78 
5.67 
0.47 
9.26 
0.47 
5.06 
0.74 
9.70 
0.95 



2.87 
0.024 
2.30 
0.024 
2.34 
0.028 
1.77 
0.017 
2.97 
0.030 



2.27 
2.27 
1.36 
1.36 
1.73 
1.73 
1.44 
1.44 
2.38 
2.38 



TABLE III. Onsite screened Coulomb integrals, Eq. (^, 
of the sp^ hybridized orbitals and s* orbital for the Si dot 
with radius R=18.9 A, InAs with R=21.1 A, and CdSe with 
R=21.1 A, in units of eV. The integrals are screened by the 
dielectric function in Eq. (|l^), which is a function of elec- 
tron-hole separation and of the radius of the quantum dot. 
For comparison, the values in parentheses are the integrals 
obtained from full screening with the dielectric constant in 
the long-distance limit e5^'(_R). 



Element {spl,spl) 



,spl) ispl,s*) 



Si 
In 

As 
Cd 

Se 



3.16 (1.26) 
1.26 (0.80) 
2.33 (1.24) 
1.86 (1.29) 
4.74 (2.51) 



1.67 (0.91) 
0.77 (0.58) 
1.33 (0.95) 
1.20 (0.99) 
2.76 (1.89) 



0.32 (0.31) 
0.24 (0.24) 
0.28 (0.24) 
0.35 (0.35) 
0.58 (0.58) 



0.26 (0.24) 
0.17 (0.14) 
0.23 (0.18) 
0.28 (0.28) 
0.48 (0.46) 



TABLE IV. Nearest-neighbor screened Coulomb integrals, 
Eq. (|l^), of the sp'^ hybridized orbitals for the Si dot with 
radius R=18.9 A, InAs with R=21.1 A, and CdSe with 
R=21.1 A, in units of eV. The integrals are screened by the 
dielectric function in Eq. (^) . The values given by the Ohno 
formula Eq. (|l^) are listed within parentheses. 





Si bonding 


Si non-bonding 


Si bonding 


2.35 (0.58) 


0.95 (0.53) 


Si non-bonding 


0.95 (0.53) 


0.55 (0.53) 




In bonding 


In non-bonding 


As bonding 


1.43 (0.54) 


0.61 (0.54) 


As non-bonding 


0.81 (0.53) 


0.42 (0.53) 




Cd bonding 


Cd non-bonding 


Se bonding 


2.41 (1.03) 


1.03 (1.02) 


Se non-bonding 


1.62 (1.02) 


0.79 (1.00) 



TABLE V. Effective masses of Si with the tight-binding 
parameters of Vogl et al. and our parameters in Table VI, 
in units of the free electron mass, rrid and rrict denote the 



longitudinal and transverse effective masses at the lowest con- 
duction energy near X. m^i and m„h are the effective masses 
at F of the two highest valence bands with a light mass and 
a heavy mass, respectively. The hole masses are averages of 
the three directions given in Ref. |38[ The cyclotron resonance 
data are taken from Ref. M 





rUci 


met 






Vogl et al. 


0.73 


1.61 


0.18 


0.39 


Our parameters 


0.91 


0.30 


0.15 


0.55 


Cyclotron resonance 


0.92 


0.19 


0.15 


0.54 



TABLE VI. Tight-binding parameters for lefectron and 
hO|l«| states of Si, in units of eV. The notationE3 of Vogl et 
alE3 is used. 







E{p) 




V{s,s) 


electron 


-3.060 


1.675 


4.756 


-8.114 


hole 


-4.777 


1.674 


8.697 


-8.465 


Vogl et al. 


-4.200 


1.715 


6.685 


-8.300 




V{x,x) 


V{x,y) 


V{s,p) 


V(s*,p) 


electron 


1.675 


21.838 


8.236 


5.994 


hole 


1.674 


4.919 


5.724 


6.133 


Vogl et al. 


1.715 


4.575 


5.729 


5.375 



TABLE VII. Tigbt-binding paraKeters for InAs, in units 
of eV. The notationtJ of Vogl et al.tJ is used. Indices a and 
c refer to anion and cation, respectively. 



E{s,a) 


E{p,a) 


E{s,c) 


E{p,c) 




-8.419 


0.096 


-2.244 


0.096 




E{s*,a) 


S(s*,c) 


V{s,s) 


V{x,x) 


V{x,y) 


12.147 


7.485 


-4.267 


1.427 


5.356 


V{sa,pc) 


V{sc,pa) 


V{s*a,pc) 


V{s*c,pa) 




4.409 


5.326 


5.846 


4.594 
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TABLE VIII. Onsite unscreened Coulomb and exchange TABLE IX. Nearest-neighbor unscreened Coulomb and 

integrals with: (0-CTO) Lowdin orthogonalized Gaus- exchange integrals with: (0-GTO) Lowdin orthogonalized 

sian-type hybrids; (NO-GTO) nonorthogonal Gaussian-type Gaussian-type hybrids; (NO-GTO) nonorthogonal Gaus- 

hybrids; and (NO-SO) nonorthogonal Slater orbitaJsi The sian-type hybrids; and (NO-SO) nonorthogonal Slater or- 

GTO integrals were calculated with the MOLPROc3 pack- bitals. rHie GTO integrals were calculated with the 

age usi«g the atomic pseudo-potentials from the Los Alamos MOLPROS p|e«kage using the pseudo-potentials from the Los 

group.tj The SO integrals are from our Monte Carlo calcula- Alamos groupCj for a two-atom molecule with a bond length 

tions. The hybrids a and b are the ones defined as sp^ and given by the bulk value. The SO integrals are from our Monte 

spl in Eq. (|lq). Carlo calculations. The indices B and A'^ designate the bond- 





O-GTO 


NO-GTO 


NO-SO 


ing and non 


bonding hybrids, respectively. 




'^Coui(a,a) of Si 


11.95 


11.65 


11.91 


Si 


0-GTO 


NO-GTO 


NO-SO 


Lu% , (a b) of Si 


9.44 


8.85 


9.00 


'^Coui(-B, B) 


8.04 


10.01 


10.60 


'^exchl",'') of Si 


1.06 


0.91 


0.73 


/ D T\T\ 


5.96 


6.65 


6.78 


Wcoui(a,a} of In 


7.90 


8.52 


7.82 


, ,0 /AT Ar\ 

'^CoullA, Aj 


4.64 


4.67 


A on 
4.89 


'^Couiia.oJ of In 


6.73 


6.54 


5.67 




0.27 


6.20 




'^cxch(",o) of In 


0.77 


0.67 


0.47 




0.11 


0.43 




/ \ r \ 

'^Coui(a,a) of As 


12.99 


12.57 


12.13 




0.04 


0.32 




^Coui(a,^') of As 


10.00 


9.54 


9.26 


T A 

InAs 


0-GTO 


NO-GTO 


NO-SO 


t^exchl",'') of As 


1.08 


0.99 


0.47 


t ID 1D\ 
'^Coull-Dl -Dj 


6.94 


8.77 


9.06 


^Coui(a>«J of Cd 


( .09 




D.09 


'^Coui(-B, A) 


5.50 


6.39 


6.59 


'^Coui(a,o) of Cd 


6.08 


5.98 


5.06 


'^Coui(A^,-B) 


5.02 


5.42 


5.43 


^c^ch{a,b) of Cd 


0.70 


0.61 


0.74 




4.12 


4.18 


4.24 


^Coui(a,a) ot be 


14.14 


lo. 1 O 


1 O Q 

iz.oo 




0.28 


4.90 




'^Coui(a,'') of Se 


10.80 


10.39 


9.70 


'^cxch(S,iV) 


0.16 


0.59 




^°xch(a,'') of Se 


1.15 


1.08 


0.90 


'^cxch(A^,A^) 


0.04 
0.04 


0.18 
0.29 






















CdSe 


0-GTO 


NO-GTO 


NO-SO 












6.94 


8.77 


9.06 










^Coui(-B, B) 


6.89 


8.66 


8.74 












5.66 


6.62 


6.84 










io°couiiN,B) 


4.85 


5.16 


5.01 












4.06 


4.11 


4.13 










^c^chiB,B) 


0.27 


4.35 












^c^ch{B,N) 


0.19 


0.69 














0.03 


0.13 












'^exch(A^,A^) 


0.04 


0.24 
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FIG. 1. Coulomb energy {eh\J{f)\eh) as a function of scal- 
ing factor / of the onsite integrals. The Coulomb energy with 
the highest hole wave function and the lowest electron wave 
function \eh} is shown for various radii of Si spherical quan- 
tum dots, with the onsite integrals scaled by the factor /. 
That is LJ ^ fu! from the values in Tables |^ and The 
offsite integrals are only indirectly scaled through the onsite 
integrals in the Ohno formula, Eq. jisl). The Coulomb energy 
is normalized by its value at / = 1. As the dot size increases, 
the Coulomb energy becomes significantly less sensitive to 
variations in the onsite integrals. 



FIG. 2. Excitonic gap of Si spherical quantum dots as a 
function of the dot radius. The photoluminescence data are 
taken from Ref. ^. The other two sets of excitonic gaps apei 
calculated with the tight-binding parameters of Vogl et alB3 
and our parameters in Table VI, respectively. Our parame- 



ters give significantly better agreement with experiment than 
the parameters of Vogl et al. This good agreement is due 
to the improved effective masses obtained with our optimized 
parameters. Inset: Coulomb shift versus the dot radius. The 
Coulomb shift does not vary much between the parameter 
sets. 



FIG. 3. Excitonic gap and single-particle gap of CdSe 
spherical quantum dots as a function of the dot radius. 
The photoluminescence excitation (PLE) gaps are taken from 
Ref. pi The scanning tunneling spectroscopy (STM) iSftpp are 
obtained from recent STM tunneling dl/dV spectraEl'S,The 
excitonic gaps of the pseudo-potential (PP) calculationalj are 
about 0.15 eV lower than the PLE gaps. Our excitonic gaps 
are in good agreement with the PLE gaps. The small differ- 
ence between our single-particle gaps and the STM quasipar- 
ticle gaps indicates that the quasiparticle polarization energy 
is small for these dots. 



FIG. 6. Unscreened exchange energy, Eq. (pl|), as a func- 
tion of cutoff distance, with the Coulomb potential replaced 
by a cutoff potential for the CdSe spherical quantum dot of 
radius _R=21.1 A. The unscreened exchange energy of four 
different types of electron-hole configurations is shown. The 
electron and hole configurations are labeled by the dominapi 
angular-momentum component of their envelope functionsO 
Except for the s-like electron and p-like hole configuration, the 
variation of the exchange interaction extends over the whole 
dot. 



FIG. 7. Unscreened exchange energy, Eq. (pl|), as a func- 
tion of cutoff distance, with the Coulomb potential replaced 
by a cutoff potential for the InAs spherical quantum dot of 
radius R=21.1 A. The unscreened exchange energy of four 
different types of electron-hole configurations is shown. The 
electron and hole configurations are labeled by the dominapt 
angular-momentum component of their envelope functionsEa 
Long-range exchange interactions appear for the s-like hole 
with both the s-like electron and the p-like electron. 



FIG. 8. Exchange charge density q(R)ch from Eq. (|2^) of 
(a) the s-like electron and s-like hole configuration, and of (b) 
the p-like electron and s-like hole configuration for the CdSe 
quantum dot with radius 21.1 A. The exchange charge density 
is plotted in a plane through the center of the dot. The unit of 
the horizontal axes is the lattice constant of CdSe. The plots 
show that the orthogonality between the electron and hole 
wave functions is global not local, with a p-like oscillation or 
a 2s-like oscillation respectively. These global oscillations of 
the exchange charge density lead to the long-range variation 
of the exchange interactions in Fig. o. 



FIG. 4. Excitonic gap and single-particle gap of InAs 
spherical quantum dots as a function of the dot radius. The 
measured PLE gaps are taken from Ref. ^. The STM-gaos are 
obtained from the tunneling spectra of Millo et aitj'Lj The 
pseudo-potential gaps (PP) are from Ref. |l^. The ap^cPa* 
tight-binding (TB) single-particl*|gaps are plotted using the 
fitting parameters of Allan et a/.tj The inclusion of d orbitals 
and spin-orbit coupling raises the gaps as much as 0.2 eV in 
comparison with our sp^ s* model. It is not understood why 
the experimental curve is so much flatter than the theoretical 
curves. 



FIG. 5. Unscreened exchange energy, Eq. (pi|), as a func- 
tion of cutoff distance, with the Coulomb potential replaced 
by a cutoff potential for various radii of Si spherical quantum 
dots. The energies are for the highest hole wave function and 
the lowest electron wave function. The curves show that there 
is an oscillation region for small cutoff distances followed by a 
saturation region beyond 15 A. This saturation suggests that 
the effective range of the exchange interaction in Si quantum 
dots is around 15 A regardless of the dot radius. 
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